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Equiatomic nickle-titanium (NiTi) is investigated to determine the consequences of point defects on the 
Martensitic phase transformation. Using molecular dynamics (MD) simulations, NiTi with 0.1%, 0.5%, 1%, 
2%, and 3% of Schottky-type defects (vacancies) have been modeled with the temperature of structural trans¬ 
formation elucidated. Increasing the concentration of point defects leads to this transformation occurring at 
lower temperatures than the perfect structure while the final monoclinic unit cell angle (7) substantially 
decreases. Modeling anti-site defects at the 0.1%, 0.5%, and 1% concentration level indicates the cubic to 
monoclinic structural transformation temperature decreases even faster with a more dramatic change in 7 
compared to the vacancy structure. The change in this Martensitic transformation stems from pinning due to 
the structural defects. 
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I. INTRODUCTION 


NiTi is an ideal material for multifarious practical 
application due to its reversible, diffusionless transfor¬ 
mation between the austenite and Martensite phases.^l 
At higher temperature, NiTi forms the cubic B2 (CsCl- 
type) crystal structure while at lower temperature the 
compound forms the monoclinic B19' crystal struc¬ 
ture, as shown in Fig-i This transformation is in¬ 
duced through external strain or changing tempera¬ 
ture and has been extensively investigated through 
computation and experiment. For example, den¬ 
sity functional (DFT) calculationPSI and molecular dy¬ 
namics (MD) simulationPI^ both accurately reproduce 
the stress-induced Martensitic transformation in NiTi. 
The temperature-induced phase transformation has also 
been analyzed by tight-binding calculations and MD 
simulation.EI On the experimental side, the temperature- 
induced transformations are often monitored via differ¬ 
ential scanning calorimetry (DSC). The shape memory 
and pseudelasticity properties of nearly equiatomic NiTi 
have been studied across a range of temperatures from 
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FIG. 1. The structural transformation of NiTi from (a) cubic B2 
structure to (b) monoclinic B19' structure with the angle, 7, 
highlighted 


-20°C to 100°C. It is shown that the loading rate as well 
as temperature could both have great effect on the shape 
memory behavior of NiTi.® 

Beyond developing a fundamental understanding of 
this structural transformation, research has also focused 
on the ability to control the transformation by vary¬ 
ing the composition. Modifying the concentration of Ni 
and Ti substantially changes the transformation temper¬ 
ature. Experiments show that increasing Ni concentra¬ 
tion leads to a decrease in Martensitic transformation 
temperature This variation in transformation tem¬ 
perature has qualitatively been confirmed using MD sim¬ 
ulations between 47% to 53% of Ni content.® 

Nearly every previous report focuses speci&ally on the 
perfect, defect-free NiTi crystal structureYet, stud¬ 
ies have shown Schottky-type defects (vacancies) as well 
as anti-site defects are both energetically favorable in 
NiTi and thus should be expected in the products .ES® 
Anti-site defects have been previously modeled through 
MD simulation to understand the crystalline to amor¬ 
phous transition of NiTi.E3 Through the random ex¬ 
change of Ni atoms and Ti atoms in the structure, the 
crystalline to amorphous transition was reported to oc¬ 
cur for a long-range order parameter (LRO) less than 
0.4. This work clearly indicates the important relation¬ 
ship between chemical disorder and phase transition. 
More recently, using a two-dimensional Lennard-Jones 
potential, the phase transformation of NiTi indicated the 
preferential formation of vacancies upon cycling. These 
can act as nucleation sites for subsequent Martensitic 
transformations .E51 

The work presented herein focuses on employing MD 
simulations to investigate the effects of point defect con¬ 
centration on the temperature induced phase transfor¬ 
mation in NiTi. Given the computational cost associated 
with DFT and the need for large supercells to capture the 
dilute concentration of defects (<3%), MD simulations 
are the most viable computational method available .E3 
Our work shows that increasing the concentration of 
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FIG. 2. The number of atoms modeled using MD for the cool¬ 
ing cycle of the NiTi Martensitic phase transformation show 
finite size effects. The bulk properties of NiTi are not reliable 
for systems smaller than 8192 atoms. 


point defects dramatically changes the temperature of 
the Martensitic phase transformation. These calculations 
highlight the importance of defect concentration on the 
defect-property relationship in NiTi. Moreover, they pro¬ 
vide a fundamental understanding of how point defects 
influence structural transformations. 


II. MODEL AND SIMULATION DETAILS 

Calculations were carried out using molecular d 3 mam- 
ics method as implemented in the LAMMPS software 
package .1^ The interatomic interactions were captured 
through the Finnis Sinclair type EAM (embedded atom 
method)!^ potential developed by Lai and Liu .23 Param¬ 
eters for this potential have been fit from first principle 
calculations of the B2 phase at 0 K. Although simpler in¬ 
teratomic potentials, e.g., Lennard-Jones, have also been 
successfully used to reproduce the phase transformation 
of NiTi,I^ EAM type potentials possess a more accurate 
description of interatomic atomic interactionsMore¬ 
over, EAM potentials account for the many body effect 
and distinguish between different coordination environ¬ 
ments, which is essential when vacancies are introduced 
into the system. For all calculations, periodic boundary 
conditions were implemented to reduce the finite size er¬ 
rors associated with the small simulation system. 

To ensure the system’s finite size does not affect the 
structural transformation, supercell calculations were 
repeated by examining various system sizes includ¬ 
ing of 1024 atoms (24A x 24A x24A), 2000 atoms 
(30 A X 30A X 30 A), 8192 atoms (48 A x 48 A x48 A), 
11664 atoms (54A x 54A x54A), and 16000 atoms 


(60 A X 60A x60 A). The transformation was determined 
by monitoring the unit cell angle (7) to indicate the for¬ 
mation of the B19' structure. As illustrated in Fig.|^ 
the transformation temperature was highest for the 1024 
atom system with a significant decrease in transforma¬ 
tion temperature observed for the 2000 atom system. 
Continuing to increase the system size to 8192 atoms 
further decreases the transformation temperature, which 
is then nearly constant with 11664 atoms. Interest¬ 
ingly, modeling 16000 atoms shows this supercell is 
large enough to accommodate the formation of nano- 
twinning domains in accordance with the experimental 
identification.ESl The transformation temperature of this 
system is still nearly the same as the 11664 atom model. 
As a result, the system with 8192 atoms was selected 
based on the compromise between accuracy and compu¬ 
tational cost. 

The initial configurations were relaxed at 400 K and 
zero external pressure for 10 ps, using a NPT ensem¬ 
ble. Nose-Hoover thermostat algorithrrPSand Parrinello- 
Rahman methods were used to control the temperature 
and pressure of the system respectively. The time step 
was set at 0.5 fs. The temperature dependent response 
of the system was examined by cooling following a ramp- 
rate of 1 K per 2000 simulation steps until a system tem¬ 
perature of 200 K was reached. A heating cycle was then 
followed using the same ramp-rate to 400 K. To inves¬ 
tigate the effect of point defects, NiTi in the B2 struc¬ 
ture was modeled by incorporating 0.1%, 0.5%, 1%, 2%, 
and 3% percent of random Schottky-type defects ensur¬ 
ing the 1:1 stoichiometry was maintained. The anti-site 
defects were also modeled as a random distribution of 
0.1%, 0.5%, and 1% disorder again making sure the el¬ 
emental ratio was constant. Visualization of the results 
was done by Ovito (Open Visualization Tool) software .E3 


ML RESULT AND DISCUSSION 
A. Structural Transformation in Perfect NiTi 

The structural relaxation of the initial, perfect NiTi in 
the austenite crystal structure (B2) yields unit cell pa¬ 
rameters at 400 K, shown in Table 1, that are in excellent 
agreement with previous reported values. ESI The relax¬ 
ation of B19' at 200 K also produces lattice param eters 
that differ by less than 3% compared to literature .^221 
The greatest differences in the optimization of B19' 
arises from variation in 7, which is likely related to the 
interatomic potential used and the boundary conditions. 
Nevertheless, the EAM used in this study is valid to de¬ 
scribe this interatomic interactions .EZI 

The structural transformation temperature is a critical 
component of shape memory alloys. The temperature 
of the phase change is often described by the starting 
(onset) temperature as well as the temperature where 
the transformation is complete, i.e., the final tempera¬ 
ture. Using this description, the transformation of per- 
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TABLE 1. The calculated NiTi lattice parameters for the B2 struc¬ 
ture at 400 K and B19' structure at 200 K. These values are 
compare d agai nst the previously reported experimentaP^ and 

calculatecP^EZl values. 




a (A) 

6(A) 

c(A) 

7(°) 

B2 


3.011 

4.258 

4.258 

90.00 

Literature 

(B2j^ 

3.013 

4.261 

4.261 

90.00 

B19' 


3.032 

4.616 

4.224 

95.30 

Literature 


2.956 

4.455 

4.189 

93.26 

Literature 


2.929 

4.686 

4.4048 

97.80 


feet NiTi is calculated to occur at Ms = 310 K, Mf 
= 305 K, As = 339 K, and Af =346K, where M is 
Martensite, A is austenite, s is the start (onset) temper¬ 
ature, and / final temperature. These calculated val¬ 
ues agree w ell w ith the reported phase transformation 
temperature This transformation also leads to the 
technologically important hysteretic behavior, which can 
be derived from cooling-heating cycle. Here, the hys¬ 
teresis is calculated by taking the difference between the 
maximum values of the derivatives for the transforma¬ 
tion (Fig.[^. In perfect NiTi, the hysteresis is 35 K for the 
8192 atom system. This is in agreement with the previ¬ 
ously calculated values, which range between 20 K and 
30 kP 


B. The Influence of Point Defects 

To understand the effects of chemical disorder on 
Martensitic transformation, the concentration of point 
defects was varied. Specifically, 0.1%, 0.5%, 1%, 2%, 
and 3% of Schottky-type vacancies were modeled as 
well as 0.1%, 0.5%, and 1% of anti-site defects. Al¬ 
though these concentrations are higher than the reported 
thermal equilibrium values,13l choosing higher concen¬ 
trations enables the effect of defects to be accentuated. 
It should be noted, however, that the synthesis proce¬ 
dure, especially when conducted at elevated tempera¬ 
tures, could yield these higher defect concentrations. 
Fig.j^demonstrates the changes in 7 with respect to sim¬ 
ulation step during cooling and heating cycles with dif¬ 
ferent concentrations of vacancies. The abrupt change of 
7 is attributed to the occurrence of the phase transforma¬ 
tion. 

A few notable features arise from the inclusion of ran¬ 
dom vacancies. Increasing their concentration leads to a 
significant change in the temperature of structural trans¬ 
formation. In all models, the increasing concentration 
of vacancies causes the transformation to occur at lower 
temperatures during the cooling cycle and the heating 
cycle. Moreover, a suppression of 7 is present at higher 
vacancy concentrations. This transformation is also cal¬ 
culated to occur over a wider temperature range indi¬ 


cated by the slope of the transformation. These changes 
are most evident in NiTi with 3% vacancies. 

The Martensitic transformations can be defined based 
on shear and dilatational displacements, the former be¬ 
ing parallel and the latter normal to the habit plane . 1 ^ 
Hence, this decrease in the temperature range and the 
supression of 7 is due to the vacancies inhibiting the 
transformation by reducing its shear component. The 
delayed phase transformation can also be explained by 
noting that the Martensite phase minimizes the internal 
energy whereas austenite maximizes the entropyP^ In¬ 
deed, introducing vacancies increases the internal energy 
and the configurational entrop}ESl favoring austenite over 
Martensite. As a result, a greater energetic driving force 
is required for the transformation to occur. 

NiTi with anti-site defect shows a similar change in 
7 versus temperature (simulation step), Fig.[^ In the 
temperature range of studied here, systems with 2 % and 
3% of anti-site defects did not show any phase trans¬ 
formation and likely occur at much lower temperature 
(<200 K) than calculated. Increasing the concentration 
of anti-site defects, the same modification of the Marten¬ 
sitic transformation occurs. The higher concentration of 
defects forces the transformation to occur at lower tem¬ 
perature and result in a smaller unit cell angle. These 
changes are even more pronounced than in the vacancy 
model with the transformation temperature with 1 % of 
anti-site defects resulting in similar properties to the 3% 
of vacancy model. Moreover, 7 in B19' is dramatically 
suppressed with an angle of ^92.5° for the 1 % anti-site 
calculation compared to 95.2° for the 1% vacancy. 

Research has shown that a slight change in the com- 
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FIG. 3. The change of unit cell angle (7) with simulation step 
and temperature (shown above) with 0%, 0.1%, 0.5%, 1%, 2%, 
and 3% of Schottky-type defects (vacancies) indicates the struc¬ 
tural transformation. The abrupt change in 7 correspond to the 
occurrence of phase change. 
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FIG. 4. The change of unit cell angle (7) with simulation step 
and temperature (shown above) with 0%, 0.1%, 0.5%, and 1% 
of anti-site defects indicates a large change in the structural 
transformation temperature. 


position of NiTi forming Ni-rich or T i-rich p hases also af¬ 
fects the phase transition drastically^EES por example, 
the addition of ^1% excess Ni causes Ms to decreases 
by lOOK.l^The average transformation temperature for 
the 1% anti-site defect structure calculated here is 81 K 
lower than the perfect structure. This phenomenon can 
be rationalized by noting that increasing the concentra¬ 
tion of Ni or Ti produces elemental islands (eith FCC or 
HCP) in the matrix. These islands energetically prefer to 
remain in the B2 (austenite) crystal structure rather than 
undergo the transformation.!^ Because the anti-site mod¬ 
els used here mimic the formation of elemental Ni and Ti 
nano-island, the suppression of transformation tempera¬ 
ture is consistent. 

The changes in the transformation temperatures with 
the incorporation of vacancy and anti-site defects are 
summarized in Fig.j^The average of the start and finish 
temperatures are plotted against the point defect concen¬ 
tration. The transformation occurs at much lower con¬ 
centrations of anti-site defects compared to vacancies. 
This analysis shows that at 1% defect concentration, the 
anti-site defects reduce the transformation temperature 
by ^26% for the cooling cycle and «19.5% for the heat¬ 
ing cycle compared to the vacancy model. Furthermore, 
the hysteresis of the transformation also decreases with 
an increasing concentration of point defects. 


IV. CONCLUSION 

Molecular dynamic simulations are a valuable tool to 
study the effect of point defects on the phase transfor¬ 
mation of NiTi during cooling and heating. Both vacancy 



FIG. 5. Ghanges of unit cell angle with simulation step and 
temperature with various percentages of vacancies (closed cir¬ 
cles) and anti-site (open squares) defects. The Martensitic and 
austenite transformations temperature are determined from 
the average of the start and finish Martensitic and austenite 
temperature respectively. 


and anti-site defects yield lower phase transformation 
temperatures and angles. The Martensitic transforma¬ 
tion changes due to pinning effects from the point de¬ 
fects. In addition, systems containing anti-site defects 
further resist to the transformation due to the forma¬ 
tion of nano-islands of the pure metals which favor the 
austenite phase. The study of the point defect within 
this phase transformation not only provides insight into 
the physics of this phase transformation but it also en¬ 
ables more fine-tune adjustment of the transformation 
temperatures through defect chemistry. 
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